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Abstract 

We present a new a-priori estimate for discrete coagulation-frag- 
mentation systems with size-dependent diffusion within a bounded, 
regular domain confined by homogeneous Neumann boundary con- 
ditions. Following from a duality argument, this a-priori estimate 
provides a global bound on the mass density and was previously 
used, for instance, in the context of reaction-diffusion equations. 

In this paper we demonstrate two lines of applications for such an 
estimate: On the one hand, it enables to simplify parts of the known 
existence theory and allows to show existence of solutions for gen- 
eralised models involving collision-induced, quadratic fragmentation 
terms for which the previous existence theory seems difficult to apply. 
On the other hand and most prominently, it proves mass conservation 
(and thus the absence of gelation) for almost all the coagulation coef- 
ficients for which mass conservation is known to hold true in the space 
homogeneous case. 

Subject Class: 35B45, 35Q72, 82D60 

Keywords: discrete coagulation-fragmentation systems, mass conservation, 
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1 Introduction 



We consider the time evolution of a physical system where a set of particles 
can aggregate into groups of two or more, called clusters, and where these 
clusters can diffuse in space with a diffusion constant which depends on their 
size. If we represent space by an open bounded set Q C with regular 
boundary, the initial-boundary problem for the concentrations q = Ci{t, x) > 
of clusters with integer size i > 1 at position x ^ Q and time t > is given 
by the discrete coagulation-fragmentation system of equations with spatial 
diffusion and homogeneous Neumann boundary conditions : 

dtCi - diA^Ci = Qi + Fi ioT X e Q,t > 0,i e W, (la) 
VxQ ■n = for X e dn,t>0,i eN\ (lb) 
Ci(0, x) = c°(x) for xeQ,ie W, (Ic) 

where n = n{x) represents a unit normal vector at a point x G dQ, di is the 
diffusion constant for clusters of size i, and 

^ i—l oo 
Qi = Qi[c] '■= Qi — Qi '■= — '^i-j ~ ^i,. 



(2) 



-^i[c] := — F- := ^ Bi^j Pi+j^i Ci+j — BiCi. 



The parameters Sj, [3ij and ajj, for integers i, j > 0, represent the total rate 
Bi of fragmentation of clusters of size i, the average number (3ij of clusters 
of size j produced due to fragmentation of a cluster of size i, and the co- 
agulation rate ajj of clusters of size i with clusters of size j. We refer to 
these parameters as the coefficients of the system of equations. They repre- 
sent rates, so they are always nonnegative; single particles do not fragment 
further, and mass should be conserved when a cluster fragments into smaller 
pieces, so one always imposes 

aM = aj,i>0, A,j>0, (z,jgN*), (3a) 
Bi = 0, Bi> 0, {i E N*), (3b) 



(ieN,i>2). (3c) 



In fact, the last condition (3c) implies the conservation of the total mass 
In Si^i ^ '^^5 which becomes obvious from the following formal fundamental 
identity or weak formulation of the coagulation and fragmentation operators: 
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Consider a sequence of nonnegative numbers {q}, and define Qj, Fi as in eqs. 
(2), then, for any sequence of numbers ipi, 



oo oo 



i=l 1=1 7=1 , , 

oo oo / i—1 \ ^ ' 

^^iFi = - Yl BiCi i^i-Yl ^i'i^i I • 

i=l 1=2 \ j=l J 

As a (still formal) consequence for solutions {q} of (1) - (2), one can calculate 
the time derivative of the integral of the moment ^ ipiCi to obtain 



^ « oo « oo 

at Jn Jn 



(5) 



since the integral of the diffusion part vanishes due to the homogeneous 
Neumann boundary condition. By choosing ipi := i above and thanks to 
(3c), we have ^^i^Qi = Yli^i'^Fi = 0, and the total mass is formally 
conserved : 

IIp(^>-)IIli = / >>Q(t,x)rfx= / >>c^(x)rfx= ||p"|L (t>0). (6) 



1= y^ici{t, x) dx = / y^ic^{x)dx 



Our main aim in this work is to provide some new bounds on the regular- 
ity of weak solutions for system (1) - (2) by means of techniques developed 
in the context of reaction-diffusion equations [9, 16, 17], and to give three 
applications to those bounds, the main one proving rigorously (for almost 
all the coefficients where this is true in the homogeneous case) mass conser- 
vation (6) and thus the absence of gelation, a well-known phenomenon in 
coagulation-fragmentation models [11, 10], where the formal conservation of 
mass is violated as clusters of infinite size are formed. 

In this paper we will work with the global weak solutions constructed in 
[15] under the assumption 

lim ^= hm ^l±l£l±i:l = o^ (for fixed z > 1), (7) 

j^+oo j j^+oo i + j 

which were later extended in [20] to the case of $7 = M^. The notion of 
solution is the following, which we take from [15]: 
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Definition 1.1. A global weak solution c = {cj}j>i to (1) - (2) is a sequence 
of functions Cj : [0, +oo) x 1] — > [0, +oo) such that for each T > 0, 



oo 



r oo 



sup 

t>0 



X) 



dx < 



X 



1=1 



dx. 



(9) 



(10) 



and for each z > 1, q is a mild solution to the i-th equation in (la), that is, 



Ci{t) = e^^^^'c'l + / e''^''^^'-'^Qi[c{s)] ds, t>0 



d,Ai{t-s). 



'11] 



where Qi[c] is defined by (2), Ai denotes the closure in L^{Q) of the un- 
bounded linear operator A of L'^{fl) defined by 

D{A) := {w e H^iCl) \Vw ■n = onOn}, Aw = Aw, (12) 

and e*^'"^^* is the Co-semigroup generated by diAi in L^{Q). 

The existence result of [15] reads: 

Theorem 1.2 (Laurengot-Mischler) . Assume hypotheses (3) and (7) on the 
coagulation and fragmentation coefficients. Assume also that 

di > for all i > 1, 

and that the non-negative initial datum has finite mass: 



cj > on f2 and 



„ oo 

Jn 



< +00. 



Then, there exists a global weak solution to the initial- boundary problem (1) 
- (2) in the sense of Definition 1.1. 

Under the extra assumptions on the diffusion constants and the initial 

data 



< inf{(ij} =: d, D := sup{di} < +oo, 

oo 



(13) 
(14) 



i=l 



we are in fact able to prove the following bound on the mass density 
p{t, x) := ^ Denoting by VLt the cylinder [0, T] x f2, we have the 
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Proposition 1.3. Assume that (3), (7), (13) and (14) hold. Then, for all 
T > the mass p of a weak solution to system (1) - (2) (given by Theorem 
1.2) lies in L'^{Qt) o,nd the following estimate holds: 

llplU^(n.) < (l + ^^^) T ||p(0, OlU^cn). (15) 

Remark 1.4. Note that the assumption (7) is only included in Proposition 1.3 
in order to ensure the existence of a weak solution via Theorem 1.2. Without 
assumption (7), the bound (15) would still hold for smooth solutions of a 
truncated version of system (1) - (2) uniformly with respect to the truncation. 
See [15] for the details of such a truncation. 

In addition to Proposition 1.3, we give a new proof of an L} bound of the 
various coagulation and fragmentation terms: 

Proposition 1.5. We still assume that (3), (7), (13) and (14) hold. Then, 
for all T > and z G N* all the terms Qf, , F^^ and associated to a 
weak solution to system (l)-(2) (given by Theorem 1.2) lie in L^{Qx) with a 
bound which depends in an explicit way on the coagulation and fragmentation 
coefficients, the diffusion coefficients, and the initial data c°. 

Remark 1.6. The fact that the terms Q^, Q~ , F^ and F[' associated to a 
weak solution are in L^{VLj') is included in the definition of weak solution; 
the main content of Proposition 1.5 is the explicit dependence of the bounds 
on the coefficients and initial data, which can be used to obtain uniform 
estimates for approximated solutions as we show for instance in section 3. 
For details on the explicit L} bounds we refer to the proof of Proposition 1.5 
in section 2. 

Remark 1.7. The bounds on , , F^ and F^ require the assumption 
(7) only to ensure existence. They would hold at the formal level (that is, for 
smooth solutions of a truncated system) under the less stringent assumption 

:= sup < +00 {le W). (16) 

isN ^ + j 

Note that the above bound also holds when assumptions (3), (7) are 
replaced by the assumptions of Theorem 1.2 in [15], but the proof is then 
much more difficult as it requires an induction on i which can be removed 
under our extra assumptions. 

In section 3, as a first application of the bounds obtained in Propositions 
1.3 and 1.5, we give a very simple proof of existence of weak solutions to 
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(l)-(2) in dimension N = 1 (that is, the result of Theorem 1.2 in dimension 
1) under the additional assumptions (3) and (7). 

Our main application of the Propositions 1.3 and 1.5 is however related to 
the problem of conservation of mass (6), which holds rigorously for solutions 
to a truncated system (see e.g [15]). Nevertheless, it is an important issue 
in coagulation-fragmentation theory whether (6) holds for weak solutions of 
system (1) - (2) itself, or if (6) is replaced by an inequahty stating that mass 
in non-increasing in time. If at some time t, the identity (6) does not hold 
any more, we say that gelation occurs, which means from a physical point of 
view that a macroscopic object has been created. 

Our main result in section 4 basically shows that (under the assumptions 
(3) and (7)) gelation does not occur when the coagulation coefficients 
at most linear and, moreover, slightly sublinear far off the diagonal i = j. 
More precisely, we prove mass conservation under the following condition on 
the coefficients aij: 

Hypothesis 1.8. There is some bounded function 9 : [0,+oo) (0, +oo) 
such that 6{x) when x —>■ +oo and 

ai,j<{i+3)e{j/i) for all j> I. (17) 

( Or equivalently, by symmetry, 

ai,j < {i+ j)9{max{j/i,i/j}) for all i,j > 1.) 

Theorem 1.9. Assume that (3), (7), (13), and (14) hold. Also, assume 
Hypothesis 1.8. Then, the weak solution to the system (1) given by Theorem 
1.2 has a superlinear moment which is bounded on bounded time intervals; 
this is, there is some increasing function C = C (T) > 0, and some increasing 
sequence of positive numbers {ipi}i>i with 

lim ipi +00 (18) 

i— >oo 

such that for all T > 0, 

1 y^ii)iCi<C{T) for allt e[0,T]. (19) 

i=i 

As a consequence, under these conditions all weak solutions given by Theorem 
1.2 of (1) conserve mass: 

I po{x) dx = pit,x) dx for all t > 0. (20) 
Jn Jn 
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Remark 1.10 (Admissible coagulation coefficients). Let us comment on Hy- 
pothesis 1.8. First note that (1.8) includes coefficients of the form 



ttij < Cst(z"/ + z^j") 

for any > such that a + jS < 1 (take 6[x) = for e > small 
enough). It is also satisfied when 

Oij < Cst 



where x i— >■ 0(x) is any positive strictly increasing function (for x big enough), 
which goes to infinity at infinity, and such that x t— > is also increasing 
(take e{X) = 0(A)-^/2)_ ^11 the examples (f) = log(l + ■), = log(l + ■) o 
log(l + -), . . . , (p = log(l + -)o- • ■olog(l + -))) satisfy this condition. Likewise, 
condition (17) also holds when (for i,j > 2) 

aij < Cst t— : h J— j — (21) 

\ \ogt log J J 

for some nondecreasing function R such that x h-> R{x)/x is nonincreasing 
and tends to when x — > +00. Note indeed that when (21) holds, 

«ij 1 R[\og{j/i) +\ogi] ^ j/i /g[logz] 



i + j logi 1 +j/ilog(j/i) +log2' 

Then, condition (17) is obtained by distinguishing the cases i > j/i and 
i j /i in both terms of the right hand side of (22). 
Assumption (21) can even be replaced by 

/.i?(log(logj)) . i^(log(logz)) 

aij <Cst\i — — — + J — — — 

V log(logz) log(logj) 

with the same requirements on R as previously. 

Note however that the linear coefficient a^j = i + j (or the coefficient 
ttij = logj + logz) does not satisfy hypothesis (1.8), though one 
would expect that Thm. 1.9 still holds for such coefficients. 



Before introducing a generalised coagulation-fragmentation model and 
thus, a third application of the Propositions 1.3 and 1.5, let us briefly re- 
view previous results on existence theory and mass conservation for the 
coagulation-fragmentation system (1). With some further restrictions on 
the coefficients as compared to [15], existence of solutions by means of L°° 
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bounds on the Cj has been proven in [3, 7, 13, 18, 19]. A different technique 
was used in [1] to prove that equation (1) is well posed, locally in time, and 
globally in time when the space dimension is one, always assuming that 
the coagulation and fragmentation coefficients are bounded. 

In a recent work [11], Hammond and Rezakhanlou considered equation 
(1) without fragmentation, and gave L°° bounds on moments of the solution 
(and as a consequence, L°° bounds on the Cj). This implies uniqueness and 
mass conservation for some coagulation coefficients that grow at most linearly 
as well as an alternative proof of the existence of L°° solutions by a-priori 
bounds on the q; for instance, if = and diffusion coefficients di are 
nonincreasing and satisfying (13) and if moreover 

oo oo 

^zc° G L-(M^), 5^^'c° G Li(M^), a,, < C [i + j) 

i=l i=l 

for some C > and all i,j > 1, then they show that mass is conserved for 
all weak solutions of eq. (1) without fragmentation. See [14, Theorems 1.3 
and 1.4] and [11, Corollary 1.1] for more details. 

In the spatially homogeneous case, mass conservation is known for general 
data with finite mass and coagulation coefficients including the critical linear 
case ttij < Cst{i + j) (see, for instance, [2, 5]). 

We finally give a third application of the Propositions 1.3 and 1.5. As 
mentioned already in the Remarks 1.4 and 1.7, Propositions 1.3 and 1.5 (de- 
spite true without restrictions on the coagulation coefficients aij for smooth 
approximating solutions) do not really improve the theory of existence of 
weak solutions for the usual models of coagulation-fragmentation like (1) as 
the full assumption (7) are needed in passing to the limit in the approximat- 
ing solutions. At best they help provide simpler proofs in particular cases, 
as done in section 3. 

On the other hand. Propositions 1.3 and 1.5 are well suited for the exis- 
tence theory of more exotic models, for instance, when fragmentation occurs 
due to binary collisions between clusters. Then, the break-up terms are 
quadratic, being proportional to the concentration of the two clusters which 
collide. This leads to coagulation-fragmentation models where all terms in 
the right hand side are quadratic. 

More precisely, we consider that clusters of size k and / collide with a 
rate bk,i > 0, leading to fragmentation. As a consequence, clusters of size 
i < maxj/c, /} are produced, in average, at a rate Pi^k,i > in such a way that 
the mass is conserved (that is, ^,j<max|fc n ^ ~ k + I). This leads to the 
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following system (for t G R+, x E Q a. bounded regular open subset of M^): 

^ oo 

dtCi - di A^Ci = 2 Cfc Q - ^ ai,k Ci Ck 

k+l=i k=l 
^ c» oo 

+ 9X1 5Z bk,iCkCipi^k,i-'^bi,kCiCk (ieW), (23) 



fc,i=l 't<max{fc,i} 



fc=l 



together with the initial and boundary conditions (lb), (Ic). For this model, 
the set of assumptions (3) is replaced by 



Pi,k,l = f3i,l,k > 0, 

bi,k = bk,i > 0, 61,1 = 0, 
iPi,k,i = k + l, 

i<max{fc,i} 



(z,jGN*), (24a) 

(z. A;,/ G < max{A;,/}), (24b) 

{i,k eN*,i < k), (24c) 

(A;,/gN*). (24d) 



Because of the quadratic character of the fragmentation terms, the in- 
ductive method for the proof of existence devised by Laurengot-Mischler [15] 
seems difficult to adapt in this case. The method presented in our first ap- 
plication can however be adapted, provided that the dimension is = 1 and 
that the following assumptions are made on the coefficients: 

Hypothesis 1.11. Assume (24), and suppose that the diffusion coefficients 
are uniformly bounded above and below (eq. (13) j and that the initial mass 
lies in L'^iVt) (eq. (14) In place of (7) we assume further that 



<^k,l 

iim — — 

>oo I 



0, 



lim sup <! (3i^k,i 

l~*oo k I Kl 



lim 

>oo I 



0, 



0. 



for fixed k eW), (25) 
for fixed z G N*), (26) 



We define a solution to (23) along the same lines as in Definition 1.1: 

Definition 1.12. A global weak solution c = {cj}j>i to (23), the bound- 
ary condition (lb) and the initial data (Ic) is a sequence of functions q : 
[0, -|-oo) X f2 ^ [0, +00) such that for each T > 0, 

c,eC{[0,T];L\Q)), t>l, (27) 

the four terms on the r.h.s. of (23) are in L^([0,T] x Q), 



sup 



„ P 00 
Jn L ,=1 



x) 



dx < 



X) 



i=l 



dx. 



(28) 



and for each z > 1, q is a mild solution to the i-th equation in (23), that is, 

c,{t) = e-^^'^i^c^ + [ e'^^^'^'-'^ Zi[c{s)] ds, t > 0, 
Jo 

where Zjc] represents the right hand side of (23) and Ai, e'^^^^^ are the same 
as in Definition 1.1. 

We are now able to prove the following theorem: 

Theorem 1.13. Under Hypothesis 1.11 on the coefficients and initial data 
of the equation, and in dimension N = 1, there exists a global weak solution 
to eq. (23) satisfying 

d G C([0,T],L^(fi)) nL2~"(l]T) {for all z G N*,T > 0,e > 0), 
for which the four terms appearing in the right hand side of (23) lie in L} (I^t) ■ 

Remark 1.14. The method of proof unfortunately does not seem to provide 
existence in dimensions N > 2. Dimension N = 2 looks in fact critical as it 
doesn't allow a-priori a bootstrap in the heat equation with right hand side 
in L^. A possible line of proof could follow [12] in the context of reaction- 
diffusion equations. In higher dimensions > 3, assuming additionally a 
detailed balance relation between coagulation and fragmentation, an entropy 
based duality method as in [9] could be used to define global weak solutions 
(see also [16]). 

Our paper is built in the following way: Section 2 is devoted to the proof 
of Propositions 1.3 and 1.5. Then Sections 3, 4, and 5 are each devoted to 
one of the three applications. In particular. Theorem 1.9 is proven in Section 
4 first in a particular case (with a very short proof), and then in complete 
generality. Theorem 1.13 is proven in Section 5. Finally, an Appendix is 
devoted to the proof of a Lemma of duality due to M. Pierre and D. Schmitt 
(cf. [17]), which is the key to Proposition 1.3. 

2 A new a priori estimate 

The solutions given in [15] are constructed by approximating the system (1)- 
(2) by a truncated system (the procedure consists in setting the coagulation 
and fragmentation coefficients to zero beyond a given finite size, and smooth- 
ing the initial data) for which very regular solutions exist. Then, uniform 
estimates for the solutions of this approximate system are proven. Finally, it 
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is shown that these solutions have a subsequence which converges to a solu- 
tion to the original system. In the proofs below it must be understood that 
the bounds are obtained for the truncated system (in a uniform way) and 
then transfered to the weak solution by a passage to the limit: the fact that 
this transfer can be done (in the case of the total mass) without replacing 
the equality by an inequality is the heart of our second application. 

We begin with the 
Proof of Proposition 1 . 3. Using the fact that 

dtp - A{Mp) = 0, inf {d,} < M{t, x) := < sup{d,}, 

we can deduce thanks to a Lemma of duality (['.), Appendix]) that p G L'^{Qt), 
and more precisely that 

l|p|U^(n.)< (l + ^^^)T||p(0,-)|U^(n), 

for all T > 0. For the sake of completeness, the Lemma is recalled with its 
proof in the Appendix (Lemma 6.2). □ 



We now turn to the 
Proof of Proposition 1 . 5. For F^' , it is clear that 

Fr <B,pe L2([0,T] xVt)(Z Li([0,T] x Vt), 

thanks to Proposition 1.3. For we use eq. (16) to write 

which is again in L^([0,T] x f2), and hence in L^([0,T] x Vt). 

For the coagulation terms, we have, since each Cj is less than p. 



1 1 ("~^ \ 



j=i \j=i 



(29) 



which is in L^([0,T] x Q) as p^ is, and the sum only has a finite number of 
terms. Finally, for we use the fact that Qf and F^ are already known 
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to be integrable: Thus, from eq. (1) integrated over [0,T] x f2, 

/ Ci{T,x)dx+ / / Q^{t,x)dxdt 
Jn Jo Jn 

< / c-(a;)rfx+ / / Ql{t,x)dxdt+ / / F^{t,x)dxdt. 
Jn Jo Jn Jo Jn 

This proves our result. □ 

3 First application: a simplified proof of ex- 
istence of solutions in dimension 1 

We begin this section with the following corollary of Proposition 1.5, in the 
particular case of dimension = 1. 

Lemma 3.1. Assume thatthe dimension N = 1, andthat (3), (13), (14) and 
(16) hold. Then, for allT>0,ie W the concentrations d e L°°([0, T] x Q) 
(where Ci are smooth solutions of a truncated version of (1) - (2), the L°° 
norm being independent of the truncation). 

Proof of Lemma 3.1. We carry out a bootstrap regularity argument. Thanks 
to Proposition 1.5, we know that (for all i G N*) 

{dt-diA)cieL\[0,T]xn). 

Using for example the results in [8], this implies that for any 6 > 0, 

ae L^-\[0,T]xn) (ieN*). (31) 

Now, eq. (31) shows that is actually more regular: from (the first in- 
equality in) (30), 

Q+ G Lt^5([0,T] X f]) for all 5 > 0, z G N*. (32) 

In addition, we already knew from eq. (29) that (for all z G N*) 

G L2([0,T] X fi), (33) 

[for which we do not need to assume that the space dimension is 1]. 
Consequently, omitting the negative terms (for alH G N*, 5 > 0), we can 
find hi such that 

{dt - diA) a < h, e d-H[0,T] X n). 
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As the Ci are positive, this imphes that 

a e LP([0, T] X n) for all p G [1, +oo[, i G N*. 
Again from (30), 

Qt G LP([0, T] X fi) for all p G [1, +oo[, i G N*. 
From this and (33), we can find hi such that 

(9i-rf,A)c, </i, GL2([0,T]xfi), 
which implies in turn that q G /^°°([0, T] x 1]) (for all i G N*). □ 

We now have the possibility to give a short proof of Theorem 1.2 in 
dimension 1 (and under the extra assumptions (13), (14)). Recall that a 
proof for any dimension can be found in [15]. 

Short proof of Theorem 1.2 in ID under the assumptions (13) and (14). 
Consider a sequence cf^ of (regular) solutions to a truncated version of sys- 
tem (1) - (2). Thanks to Proposition 3.1, we know that for each i G N*, 
supj^^ ll'^i^^llL°°(f7T) ^ +00. Then (for each i G N*) there is a subsequence 
of the (c*^)mgn (which we still denote by (cf^)MGN), and a function q G 
L°°(1]t), such that 

cf ^ Ci weak-* in L°^{VLt)- (34) 

Using Proposition 1.5, we also see that (for any fixed i G N*), the L}{VLt) 
norms of C^'^\ C^'^ , F^'^^ , F['^^ (the coagulation and fragmentation terms 
associated to {cf^}) are bounded independently of M. Using eq. (la) and the 
properties of the heat equation, one sees that for each i G N*, the sequence 
{cf^} lies in a strongly compact subset of L^{VLt)- Hence, by renaming our 
subsequence again, we may assume that 

cf ^ d in L^{VLt) strong , for all z G N*. (35) 

In order to prove that {cj} is indeed a solution to eq. (1) - (2), let us 
prove that all terms F^'^'\ Fr'^\ C+'^^ Cp^^ converge to the corresponding 
expressions for Cj, which we denote by Fp, F~ , , C~ , as usual. 

1. Positive fragmentation term: for each fixed i, the sum 

— 2^ ^i+j Pi+j,i ^i+j 

13 



converges to in L}{VLt) because the tails of the sum converge to 
uniformly in M (this is due to hypothesis (7)): 



'0 Jn 




dxdt <2 ( sup 

J>Jo 



Bi+j f^i+j,i 



i + j 



+ sup \\Ci+j -Ci+jllmnr)- 
j<Jo 



2. The negative fragmentation term is just a multiple of cf , so the con- 
vergence in L^{Qt) is given by (35). 

3. For each fixed i, the positive coagulation term is a finite sum of terms 
of the form aijc^^cj'' . Thanks to (34) and (35), this converges to aijCiCj 

4. The negative coagulation term is 



Since cf^ converges to q weak-* in L°°{Qt), it is enough to prove that 
YlJLi C'ij converges to YlJLi C'ij strongly in L^{Qt)- Observing 
that 



'0 Jn 




dxdt < 2 ( sup 

J>Jo 



p+SUp ||cf -Cj||ii(nj,), 
j<Jo 



we see thanks to (7) and (35) that this convergence indeed holds. 



□ 



4 Second application: mass conservation 

We begin this section with a very short proof of Theorem 1.9 in a particular 
case in order to show how estimate (15) works. More precisely, we consider 
the pure coagulation case with aij = ^/Tj and Bi = (no fragmentation), 
and with initial data satisfy additionally X^^g ^ ^^S^ Cj(0, x) dx < +oo (which 
is sightly more stringent than only assuming finite initial mass). 
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Then, using the weak formulation (4) with = log(z) (and remembering 
that log(l + x) < Cst ^/x) 

— / ^ilogicidx = / Ci Cj ( i log(l + -) + j log(l + -)]dx 

Jn -^^ Jn -^-^ \ ^ J / 

„ oo oo „ 

<2 / ^^ijCiCjdx <2 / p{t,xfdx. (36) 
Jn ^.^1 Jn 

As a consequence, we have for all T > 

« oo „ oo „7^ „ 

/ z log z Cj(T, x) (ix < / z logz Cj(0, x) (ix + 2 / / pit^xYdxdt, 

which ensures the propagation of the moment /^^q« logiCi{-,x)dx, and 
therefore gives a rigorous proof of conservation of the mass for weak solutions 
of the system: no gelation occurs. 

Our general result is obtained through a refinement of this argument 
under hypothesis (1-8). Before giving the proof of Theorem 1.9 we need two 
technical lemmas, which will substitute the intermediate step in (36). 

Lemma 4.1. Let {fii}i>i and {z/j}j>i be sequences of positive numbers such 
that {yUj} is bounded, 



Ejii = +00 and lim z/j = +oo. 
j— >+oo 

i=l 

Then we can find a sequence {C,i}i>i of nonnegative numbers such that 

oo 



-cxo, 

i=l 

i 



Ci ^ /^t ^''T'd ipi := — f^"^ all i > 1. 

Proof. We may assume that Ui is nondecreasing, for otherwise we can consider 
Ui := infj>j{z/j} instead of Ui. Then, in order to find ^ it is enough to define 
recursively := and, for i > 1, 

jji if /ii + < i'i) 

otherwise. 
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By construction, C,i < /^i for alH > 1, and also Ylj=i — ^'^^ i > 1, as we 
are assuming {z/,} nondecreasing. 

To see that {^i} cannot be summable, suppose otherwise that Xli^i 6 = 
S < +00. Take a bound M > of {/ij}, and choose an integer k such that 
Ui > S + M for all i > k. Then, by definition, 

= fii for all i > k, 

which implies that {^j} is not summable, as {/ij} is not, and gives a contra- 
diction. □ 

Lemma 4.2. Assume (17). There is a nondecreasing sequence of positive 
numbers {ipi}i>i such that ipi — > +oo when i — > +cx), and 

ttiji^i+j - ^pi) < Cj for all i, J > 1, (37) 

for some constant C > 0. 

In addition, for a given sequence of positive numbers Xi with limj^+oo K = 
+00, we can choose ipi so that ipi < Aj for all i. 

Proof. First, we may assume that the function 6 given in Hypothesis 1.8 is 
nonincreasing on [1, +00), as we can always take 9{x) := sup^^^ 9{y) instead. 

We choose a sequence of nonnegative numbers {^j} by applying Lemma 
4.1 with 




min < Ai, — -=-, } . (39) 



Note that the conditions in Lemma 4.1 are met: the sequence in the right 
hand side of (38) is not summable, and the right hand side of (39) goes to 
+00 with i. If we define ipi := ^2]=! 0' then the following is given by Lemma 
4.1: 

< , , .^ , V'i < ^, A <K i> 1, 

(l + z)log(l + z) ^(\A/2) 



lim ipi = +00. 

j— >+oo 

These conditions essentially say that ipi grows slowlier than loglog(z), slowlier 
than 6{^^/i/2)~^ , and slowlier than Aj, yet still diverges as i — > +00. 

We can now prove (37) to hold for these {ipi} by distinguishing three 
cases: 
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1. For any i,j > 1, as log(l + k) > 1/2 for all A; > 1 

i+j i+j 



^,+,-^, = J2 ^k<2j2 ^<21og(^ + j + l)-21og(z + l)<^. 

k=i+l k=i+l 



Then, in case j < i we use the fact that 6{x) < Cq for some constant > 
and all a; > and have 

ai,j{^i+j - ^i) < 2Ce(i + j)- < 4Ce j, for j < i. 



,•2 

1 



2. Secondly, ioi i < j < i 

< E < E (A: + l)log(A:+l) 

< loglog(2^2 + 1) - loglog(2 + 1) < log (" ^Mv^^A < 

Viog(z + 1) y 

for some number Ci > 0. Thus, 

aijiipi+j - ^i) < CiCeii+j) < 2CiCej. 

3. Finally, for j > i"^, 

and as 6 is nonincreasing on [1, +oo) (we may assume this; see the beginning 
of this proof), we have for all j > 

Together, these three cases show (37) for all i, j > 1. □ 

Now we are ready to finish the proof of our result on mass conservation: 

Proof of Theorem 1.9. As remarked above (cf. beginning of section 2), we 
will prove the estimate (19) for a regular solution to an approximating system, 
with a constant C{T) that does not depend on the regularisation. Then, 
passing to the limit, the result is true for a weak solution thus constructed. 

We consider a solution to an approximating system on [0, +oo), which 
we still denote by {ci}i>i. Then, by a version of the de la Vallee-Poussin's 
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Lemma, (see, for instance. Proposition 9.1.1 in [1] or also proof of Lemma 
7 in [6]), there exists a nondecreasing sequence of positive numbers {\i}i>i 
(independent of the regularisation of the initial data) which diverges as i — > 
+00, and such that 



/ y ^ i AiC° dx < +00. 
i=i 



(40) 



If we define rj := idj, note that this is just the claim that one can find Aj 
as above with Ajr^ < +oo. 

Taking {ipi} as given by Lemma 4.2, such that ipi < Aj for all i > 1, we 
have thus Yli^i i ipi c^i{.x) dx < +(X). Then, as integrating over VL makes the 
diffusion term vanish due to the no-flux boundary conditions, we estimate 



d 
It 



/ y^iiiiCidx < - / aijCiCjd; 



i+j)'4'i+j-i'>Pi-j'^j)dx, (41) 



where we used that the contribution of the fragmentation term is nonpositive, 
as can be seen from (4) with ifi = iipi, and the fact that 

i-l i-1 

^ j'^j < X] = ^ 

as Ipi is nondecreasing and (3c) holds. Continuing from (41), by the symmetry 
of the Oij , and using the inequahty (37) from Lemma 4.2, we have 



d 
dt 



« oo „ oo „ 

/ iipiCidx < I ai j Ci Cji {ipi+j — ipi) dx < C / dx. (42) 



Thus, Proposition 1.3 showing p G L'^{Qt) proves that ^ Q da; is 

bounded on bounded time intervals. Mass conservation is a direct conse- 
quence of this. □ 

Remark 4.3 (Absence of gelation via tightness). It is interesting to sketch an 
alternative proof showing conservation of mass via a tightness argument and 
without establishing superlinear moments. By introducing the superlinear 
test sequence i(pk{i) with (pk{i) = ^i^i<k + h>k for all k G N*, we use 
the weak formulation (4) to see (as above) that the fragmentation part is 
nonnegative for superlinear test sequences, and use the symmetry of the Qij 
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to reduce summation over the indices i > j eW, which leads to the estimate 
j Y^c, i<Pk{i) dx< f ^^'J- ['^^] [^j] ( ^^<^ 

i=l "'^ i>j j=l \ 6V / 



J 



log(l + ^) ^ log(f 



For the first term, we use log(l + j/i) < j/i- Then, for the second and third 
terms, we distinguish further the areas where i/j < log{k) and i/j > \og{k). 
When i/j < \og{k), we estimate 1 + i/j = 1 + i/j < 1 + log(/c) and k/j < 
1 + i/j < 1 + log(fc), respectively. On the other hand, when i/j > log(/c), 
both the second and the third term are bounded by one. Altogether, we get 
thanks to assumption (1.8), i.e. ^ < Cst 0{i/ j) for i < j: 

1*^-'] [^'^J'] ~r IVi>log(fc);i<fc 



iOgl«;j i/j>log(fc) V J , 




and the right hand side tends to zero a.s k oo. Hence, using Proposition 1.3 
and integrating over a time interval [0,T], we get thanks to a tightness ar- 
gument that the mass is indeed conserved, and no gelation occurs. 



5 Third Application: Fragmentation due to 
collisions in dimension 1 

Proof of Theorem 1.13. We introduce {cf^)M a sequence of smooth solutions 
for a truncated version of eq. (23). We first observe that Proposition 1.3 
still holds thanks to the duality estimate, that is p := Xlj^Cj G L'^iQ.T) for 
all T > 0. Estimate (30), in which only the coagulation kernel appears, also 
holds. Moreover, thanks to (24d), 

&fc,/CfcQ/3ifc; < Csti^^(A; + /)CfcQ < Cstip^ G L^(fiT)- 
k,l max{k,l}>i k I 
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The loss terms 



Ck, 



k=l 



oo 
A;=l 



lie then in L^[Qx) by integration of the equation on [0,T] x Q. 

Using now eq. (23), we see that (for all i G N*) dtcf — did^xcf belongs to 
a bounded subset of L^{Qt)- As a consequence, cf^ belongs (for all i G N*) 
to a compact subset of L^~^([0,T] x Q) for all T > and e > 0. We denote 
(for all i G N*) by q a limit (in L^~^([0,T] x f2) strong) of a subsequence of 



)M£N 



(still denoted by (c: 



M/eNj 



We now pass to the limit in all terms of the r.h.s. of eq. (23). The first 
term can easily be dealt with, since it consists of a finite sum. Then, we pass 
to the limit in the second term: 




Jn 



< 



'0 Jn 




K 



^ ^ (^i,k Q Cfc 
k=l k=l 
K 



dxdt 



k=l 



k=l 



dxdt 



L2 sup 

k>K 



The second part of this expression is small when K is large enough thanks 
to assumption (25), (26), while the first part tends to for all given K. 

The fourth term of the r.h.s. of eq. (23) can be treated exactly in the 
same way. We now turn to the third term: 




< 



Jn 
'0 Jn 




E E ^'^'^ ^'''''^ ~ E E ^'''^ ^'''^'^ 

fc,i=l i<max{fc,i} k,l=l 't<max{fc,i} 

ft: k<K,l<K K k<K,l<K 

E E ^'^'^ "^fc ^^''^'^ ~ E E ^'^■'^ ^'^ 
A:,/=l j<max{fe,/} k,l=l j<max{A;,/} 

bk,i 



dxdt 



dxdt 



+ 4 



^2 sup sup 



Once again, the second term is small when K is large enough thanks to 
assumption (25), (26), while the first term tends to for all given K. □ 
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6 Appendix: A duality lemma 

We recall here results from e.g. [17, 9]. We start with the 
Lemma 6.1. Assume that z : [0, +oo) satisfies 

dtz + MAz = -H on Q, 



\/z-n = on dn, (43) 

z{T, x) = on Q, 



where H G L'^inr), and di > M > do > 0. Then, 



N(o,-)llL^(n)< + (44) 



Proof of Lemma 6.1. Calculating the time derivative of | V;z|^, or alterna- 
tively multiplying eq. (43) by Az and integrating on we obtain 



1 d 
'2~dt 



I \Vzfdx+ I M{Azfdx= f -HAzdx, 
Jn Jn Jn 



where the boundary condition on z was used. Integrating on [0, T] and taking 
into account that z(T,x) = 0, 

- / \Vz{0,-)f dx + I M{Azfdxdt= f HAzdxdt 

<\\H\\mn^)\\^4mnT)- (45) 
Using that M > do we see that M(Az)^ > do ll^^llLam^), so (45) implies 
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From this and (43) we have 



Finally, 



□ 



Lemma 6.2. Assume that p : Qt [0, +oo) and satisfies 

dtp - A{Mp) < onn, (46) 
V{pM) ■n = on on, 

where M : Qt is a function which satisfies di > M > do > for some 
numbers di, do. Then, 

IIpIIl^(^.)< (i + l) TMo,-)h. 

Proof of Lemma 6.2. Consider the dual problem (43) - (44) for an arbitrary 
function H G L'^{VLt), with H > 0. Then, z > 0, and integrating by parts in 
eq. (43), one finds that 



pH dxdt = - p{dtz + MAz) dxdt 

z{dtp - A{pM)) dxdt + / p{0,-) z{0,-)dxdt < / p{0, ■) z{0, ■) dxdt, 

where we have used eq. (46), eq. (44) and the boundary conditions on pM 
and z. Hence, for any nonnegative function H G L'^{Qt), 



IT 

and thanks to Lemma 6.1 



[ pHdxdt<\\piO,-)\\^,^^^\\ziOr)\\ 



L2(n) ' 



/ pHdxdt < (1 + di/do) T ||p(0, ■) ||^.(^) ||i/|L2(n,) • 

Remembering that p > 0, we obtain by duality: 

||plL2(^^)<(l + rfi/rfo)r||p(0,-)IL2(n). 
This proves the lemma. □ 



22 



References 

[1] H. Amann. Coagulation-fragmentation processes. Archive for Rational 
Mechanics and Analysis, 151(4):339-366, 2000. 

[2] J. M. Ball and J. Carr. The discrete coagulation-fragmentation equa- 
tions: Existence, uniqueness, and density conservation. Journal of Sta- 
tistical Physics, 61(l):203-234, 1990. 

[3] P. Benilan and D. Wrzosek. On an infinite system of reaction- 
diffusion equations. Advances in Mathematical Sciences and Applica- 
tions, 7(l):351-366, 1997. 

[4] J. A. Caiiizo. Some problems related to the study of interaction ker- 
nels: coagulation, fragmentation and diffusion in kinetic and quantum 
equations. PhD thesis, Universidad de Granada, June 2006. 

[5] J. Carr and F. Costa. Asymptotic behavior of solutions to the 
coagulation-fragmentation equations. II. Weak fragmentation. Journal 
of Statistical Physics, 77(1):89-123, 1994. 

[6] J. Carrillo, L. Desvillettes, and K. Fellner. Exponential decay towards 
equilibrium for the inhomogeneous Aizenman-Bak model. Communica- 
tions in Mathematical Physics, 278(2) :433-451, 2008. 

[7] J. F. Collet and F. Poupaud. Existence of solutions to coagulation- 
fragmentation systems with diffusion. Transport Theory and Statistical 
Physics, 25(3):503-513, 1996. 

[8] L. Desvillettes and K. Fellner. Entropy methods for react ion- diffusion 
equations: slowly growing a-priori bounds. Revista Matemdtica Ibero- 
americana, 24(2):407-431, 2008. 

[9] L. Desvillettes, K. Fellner, M. Pierre, and J. Vovelle. About global ex- 
istence for quadratic systems of reaction-diffusion. Journal of Advanced 
Nonlinear Studies, 7(3):491-511, 2007. 

[10] M. Escobedo, P. Laurengot, S. Mischler, and B. Perthame. Gelation 
and mass conservation in coagulation-fragmentation models. Journal of 
Differential Equations, 195(1):143-174, 2003. 

[11] M. Escobedo, S. Mischler, and B. Perthame. Gelation in coagulation 
and fragmentation models. Communications in Mathematical Physics, 
231:157-188, 2002. 



23 



[12] T. Goudon and A. Vasseur. Regularity analysis for systems of reaction- 
diffusion equations, to appear in Annates de I'Ecole Normale Superieure. 

[13] F. Guia§. Convergence properties of a stochastic model for coagulation- 
fragmentation processes with diffusion. Stochastic Analysis and Appli- 
cations, 19(2):245-278, 2001. 

[14] A. Hammond and F. Rezakhanlou. Moment bounds for the smolu- 
chowski equation and their consequences. Communications in Mathe- 
matical Physics, 276(3) :645-670, 2007. 

[15] P. Laurengot and S. Mischler. Global existence for the discrete dif- 
fusive coagulation- fragmentation equation in L} . Revista Matemdtica 
Iberoamericana, 18:731-745, 2002. 

[16] M. Pierre. Weak solutions and supersolutions in for reaction-diffusion 
systems. Journal of Evolution Equations, 3(1):153-168, 2003. 

[17] M. Pierre and D. Schmitt. Blowup in reaction-diffusion systems with 
dissipation of mass. SIAM Review, 42:93-106, 2000. 

[18] D. Wrzosek. Existence of solutions for the discrete coagulation- 
fragmentation model with diffusion. Topological Methods in Nonlinear 
Analysis, 9(2):279-296, 1997. 

[19] D. Wrzosek. Mass- conserving solutions to the discrete coagulation-frag- 
mentation model with diffusion. Nonlinear Analysis, 49(3):297-314, 
2002. 

[20] D. Wrzosek. Weak solutions to the cauchy problem for the diffusive 
discrete coagulation-fragmentation system. Journal of Mathematical 
Analysis and Applications, 289(2):405-418, 2004. 



24 



